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SUMMARY 


Studies  in  the  Aerospace  Medical  Research  Laboratory  utilize  radioactive 
tracer  microspheres  to  determine  blood  flow  to  a specific  region,  tissue  or  organ 
under  +GZ  acceleration.  The  microspheres  are  injected  and  become  distributed 
throughout  the  body  similar  to  red  blood  cells.  They  become  trapped  in  capillary 
beds  according  to  size.  Tissue  samples  are  taken  and  radiation  measured  in  a 
gamma  counter. 

This  report  describes  a computer  program  developed  to  calculate  corrected 
counts  per  minute,  regional  blood  flow,  and  cardiac  output  from  raw  data.  The 
program  offers  a rapid,  easily  used  method  for  processing  data  obtained  from  the 
gamma  counter. 
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INTRODUCTION 


Studies  similar  to  those  of  Sostre  (7)  utilizing  plastic  radioactive  tracer 
microspheres  to  calculate  blood  flow  to  a specific  region,  tissue,  or  organ  under 
+GZ  acceleration  are  being  done  in  this  Laboratory.  The  microspheres  are  injected 
into  the  left  ventricle,  and  distributed  throughout  the  body  in  a manner  similar 
to  red  blood  cells.  They  are  entrapped  in  specific  capillary  beds  according  to 
blood  distribution  and  the  size  microsphere  used.  Tissue  samples  from  areas  of 
Interest  are  weighed  and  placed  in  a gamma  counter  for  determination  of  microsphere 
density  in  terms  of  counts  per  minute  (CPM) . A computer  program  was  developed  to 
calculate  corrected  CPM,  regional  blood  flow,  and  cardiac  output  from  raw  data. 

This  program  offers  a rapid,  easily  used  method  for  processing  data  obtained  from 
the  gamma  counter. 

Section  2 is  a brief  outline  of  the  experimental  design.  The  theory  behind 
the  program  is  presented  in  Section  3,  and  Section  4 contains  the  execution  of 
the  program.  Section  5 includes  the  program  listing. 


EXPERIMENTAL  PROCEDURE 

The  animals  used.  In  this  Instance  baboons,  are  lightly  anesthetized  with 
ketamine.  Two  femoral  artery  cutdowns  are  made;  a catheter  is  inserted  into  the 
left  ventricle  via  one  femoral  artery.  Microsphere  injections  are  made  through 
this  catheter  using  a remotely  controlled  microsphere  injection  system  (2) . A 
small  silicone  coated  polyethylene  catheter  is  inserted  into  the  other  femoral 
artery.  It  is  placed  with  the  tip  in  the  descending  aorta  just  below  the  arch. 
Fifteen  seconds  prior  to  the  injection  of  microspheres,  a constant  rate  reference 
sample  is  withdrawn.  A Harvard  Pump  @ (7)  is  used  in  this  Laboratory,  but  any 
constant  rate  withdrawal  system  is  acceptable.  (6)  Withdrawal  of  blood  continues 
throughout  the  injection  series  until  one  minute  after  the  final  injection.  The 
volume  of  blood  withdrawn  is  calculated,  and  radioactivity  is  determined.  Blood 
flow  to  this  reference  sample  is  taken  to  be  that  of  a "known"  organ;  flow  to  any 
other  organ  then  be  calculated  by  relating  to  the  blood  flow  of  the  reference 
sample  (1) . 

The  animal  is  placed  in  the  Oloff  Primate  Restraint  System  (PRS)  (5),  and 

allowed  to  recover  from  the  ketamine.  The  PRS  is  then  placed  on  the  Dynamic 

Environment  Simulator  (DES),  a large  man-rated  centrifuge  located  in  this 

Laboratory.  When  the  animal  is  fully  awake,  the  acceleration  profile  begins. 

Peak  acceleration  is  either  at  +3GZ  or  +5GZ,  with  an  equal  number  of  animals 

exposed  to  each.  Five  types  of  15  micron  radiolabeled  microspheres  are  used  hr 

each. animal:  Cerium  141,  Chromium  51,  Strontium  85,  Scandium  46,  and  Iodine  125. 

0.025  to  0.050  microcuries  are  injected,  depending  on  the  relative  activity  of 

the  isotope.  The  first  isotope  is  injected  at  +1GZ  and  is  used  to  establish  a 

baseline.  When  peak  acceleration  is  reached,  the  remaining  isotopes  are  injected 

at  30  second  intervals.  After  injection  of  the  last  isotope,  peak  acceleration 

is  maintained  for  30  seconds  bringing  total  -KS  exposure  time  to  120  seconds. 

z 

When  the  acceleration  profile  is  completed,  the  animal  is  sacrificed  with 
Uthanol(R}  and  the  brain  is  removed  and  fixed  in  formalin  for  72  hours.  One  half 
of  the  brain  is  sectioned  into  various  regions,  such  as  brainstem,  cerebellum, 
and  hypothalamus.  The  other  half  is  placed  in  several  vials  and  used  to  estimate 
total  cerebral  blood  flow.  Tissue  samples  are  also  taken  from  various  organs 
throughout  the  body,  with  each  location  assigned  a specific  ID  number.  All 
samples  are  weighed  and  placed  in  an  appropriate  vial  to  be  counted  in  a Packard 
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Auto-Gamma  Scintillation  Spectrometer Ql  When  a sample  is  too  large  to  fit  into 
one  counting  vial,  it  can  be  sectioned  and  placed  into  several  vials.  All  the 
vials  can  then  be  incorporated  into  a group  and  treated  as  one  sample.  Information 
from  the  gamma  counter  is  printed  out  on  a standard  teletype  with  a punch  tape 
capacity. 

Disposal  of  radioactive  animal  carcasses  is  in  accordance  with  Appendix  B, 
Title  10,  Code  of  Federal  Regulations,  Part  20. 


THEORY 


Program  GAMMA  was  written  for  a Control  Data  Corporation  6600  computer 
located  at  the  Aeronautical  Systems  Division  Computer  Center  at  Wright-Patterson 
Air  Force  Base,  Ohio. 

There  are  three  steps  involved  in  using  this  program: 

1.  Create  a permanent  file  of  data  from  the  gamma  counter.  This  file 
will  contain  counts  of  all  reference  samples,  standards,  assay 
samples,  and  residue  samples.  Standards  are  aliquots  of  the  micro- 
sphere solution  of  the  same  age  as  those  injected  into  the  animals. 
Residue  samples  include  samples  from  the  withdrawn  catheters,  post- 
experiment syringes,  injection  catheters,  and  manifold  washouts. 

2.  Run  program  GAMMA  using  the  above  file,  which  is  given  the  local 
file  name  Tape  1.  A local  file  named  Tape  2 will  be  created, 
this  file  contains  the  corrected  CPM,  regional  blood  flow,  and 
cardiac  output. 

3.  Catalog  Tape  2. 

Computation  of  corrected  CPM,  regional  blood  flow,  and  cardiac  output  is  as 
follows: 

1.  Counts  from  all  samples  are  corrected  for  background  radiation: 


where  c^  = original  count,  window  i 

e.^  “ mean  count  window  i,  from  empty  vials 

c'j"  background  corrected  count,  window  i 

2.  Counts  made  on  days  previous  to  the  assay  counts  are  corrected  for 
radioactive  decay  during  the  Intervening  time: 

decay  corrected  count  - original  count  x 

where  c ■ ln(.5) 

t ■ days  elapsed  since  count 
h - half-life  of  given  isotope,  in  days 


life  values  (4) 

used  are: 

Iodine 

60  days 

Cerium 

32  days 

Chromium 

28  days 

Strontium 

64  days 

Scandium 

84  days 

Counts  of  standards  are  corrected  for  decay  using  the  same  half-life 
for  each  window;  the  isotope  is  assumed  to  decay  at  the  same  rate  over 
the  entire  spectrum.  Counts  of  samples  that  are  not  standards  are 
corrected  for  decay  using  the  half-life  of  the  corresponding  isotope 
for  each  window;  the  decay  rate  within  each  window  is  assumed  to  be 
the  decay  rate  of  the  predominant  isotope. 

3.  The  total  amount  of  each  Isotope  is  computed  as  follows: 


= yi  x volume  injected  (ml) 
volume  counted  (ml) 


where  *=  total  amount  of  isotope  i injected,  not  corrected  for 

residue  samples 

n. 


± £ 

ni 

1 k-1 


5 

£ 

j"l 


'ijk 


where  y ■ average  count  of  the  standard  for  isotope  i 

n^  = the  number  of  times  the  count  for  isotope  i is  repeated 

c^jk  * the  count  of  the  standard  for  isotope  i in  the  jth  window 
repeated  the  kth  time,  corrected  for  decay  and  background 
radiation 


Thus  y^  is  the  average  of  the  total  counts  In  all  windows  from  isotope  i. 

4.  Corrected  counts  (corrected  for  overlapping  windows)  are  computed  by 
the  formula*: 

T - (F'F)”1  F'C 

where  T is  a 5x1  matrix;  t^  - corrected  counts  from  Isotope  1 

F is  a 5x5  matrix;  f ^ » proportion  of  isotope  i in  window  j 

* Rosenblatt,  Judah.  Personal  communication 
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C is  a 5x1  matrix;  cj  * observed  count  in  window  j,  corrected 

for  decay  and  background  radiation 

F is  estimated  from  the  counts  of  the  standards  as  follows: 


i V- 

ni  k-1 


i ni 


£ n £ 
i=l  k-1 


where  n^  and  c^^  are  as  in  3 above 

5 . Regional  blood  flow  is  computed  by  the  formula  (5) : 

t x withdrawal  rate  (ml/min) 

H " 

r^  x weight  of  sample  (grams) 

where  g^  * regional  blood  flow  imnedlately  after  Injection  of  isotope  i,  in 
[ml/min/gram] 

t^  - the  corrected  count  from  isotope  i 

r^  - the  sum  of  corrected  counts  from  isotope  i,  summed  over  all 
reference  samples 

6.  Cardiac  output  Immediately  after  inlection  of  Isotope  i is  given  by  the 

, . w.  x withdrawal  rate  (ml/min) 

formula  (3):  s.^  - i 


where  w^  - the  total  corrected  count  of  isotope  i injected,  computed  by 
subtracting  corrected  counts  of  all  residue  samples  from  the 
total  count  of  isotope  1 (z^  in  3 above) . 


EXECUTION  OF  PROGRAM  GAMMA 


The  permanent  data  file  described  In  the  theory  section  should  contain  the 
following  Information  for  each  sample: 

Line  1 Columns  1-6:  sample  Identification  number 

Lines  2-6  Counts  for  each  of  the  five  Isotopes  in  ascending  order 
by  energy  level 

Columns  15-20 : energy  lev*.’ , center  of  window 


Columns  50-55:  count  for  given  window 

The  sample  identification  number  should  be  unique;  no  two  samples  should  have  the 
same  ID,  with  the  exception  of  multiple  counts  of  the  standards.  The  energy  level 
at  the  center  of  the  window  should  be  within  the  following  ranges: 


Iodine 

Cerium 

Chromium 

Strontium 

Scandium 


29-35 

142-148 

317-323 

511-517 

1117-1123 


The  above  format  is  the  standard  output  of  the  Packard  Auto-Gamma  Scintillation 
Spectrometer  R (see  Appendix  A) . The  program  can  be  changed  to  accomodate  the 
use  of  a counter  with  a different  format. 

If  the  data  from  an  experiment  are  on  several  different  paper  tapes  they  should  be 
read  into  the  computer  and  combined  into  one  permanent  file  using  Editor. 

Execution  of  program  GAMMA  is  as  follows: 

1.  Attach  the  permanent  file  of  data  from  the  gamma  counter 
ATTACH,  TAPE1,  PFN ,CY-XX, ID-XXXXXX 

2.  Request  permanent  file  space  for  output  file 
REQUEST,  TAPE2,  *PF 

3.  Attach  the  permanent  file  which  contains  program  GAMMA 
ATTACH,  GAMMA,  WARD35.CY-1.ID-L740530 

4.  Attach  the  subroutine  library  IMSL* 

ATTACH , IMSL , ID-LIBRARY , SN-ASD , MR-1 


* International  Mathematical  & Statistical  Libraries.  ASD  Computer  Subprogram 
Library 
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5.  Make  IMSL  part  of  the  local  library 
LIBRARY,  IMSL 

6.  Compile  program  GAMMA 
FTN(I-GAMMA,L-0) 

7 . Execute  the  program 
LGO 

At  this  point,  the  program  will  ask  for  certain  information  which  should  be  typed 
in  as  requested.  Appendix  B contains  an  abbreviated  sample  with  typed-in 
information  underlined. 

8.  Catalog  Tape  2 

CATALOG,  TAPE2,PFN,CY-XX,RP-999 

where  PFN  is  the  permanent  file  name  and  XX  is  the  cycle  number.  A sample  of  the 
output  from  Tape  2 is  included  in  Appendix  C. 

NOTES: 

1.  Entering  sample  ID  numbers 

a.  A group  of  numbers  may  be  entered  by  putting  a minus  sign  in  front  of  the 
selected  number  satisfying  the  range.  For  example,  "580  -592"  specifies 
all  thirteen  ID  numbers  between  and  including  580  and  592.  It  is  important 
to  leave  a space  between  the  first  number  and  the  minus  sign. 

b.  If  too  many  numbers  are  specified,  the  extra  numbers  will  be  dropped 
starting  at  the  end  of  the  list. 

c.  If  too  few  numbers  are  specified,  the  program  will  not  continue  until  more 
numbers  are  given. 

4.  Do  not  include  ID  numbers  in  the  list  which  do  not  exist  in  the  file 

Taps  1.  For  example,  suppose  a group  consists  of  samples  33  -45,  and  47 
and  there  is  no  sample  46  in  the  data  file.  If  the  ID  numbers  for  the 
group  are  specified  as  "33  -47"  instead  of  "33  -45  47"  sample  47  will  be 
dropped  from  the  group  because  too  many  ID  numbers  are  specified. 

2.  Samples  counted  previous  to  assay  are  grouped  by  the  number  of  days  elapsed 
from  the  count  of  the  given  sample  to  the  assay  counts.  Any  number  of  groups 
may  be  specified.  After  all  groups  have  been  entered,  enter  0 for  the  number 
of  samples,  and  0 for  the  number  of  days  elapsed. 
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When  sample  ID  and  weight  are  requested  for  each  sample,  the  sample  ID  to 
enter  is  the  identification  of  the  location  in  the  animal  from  which  the  sample 
came. 

Assay  samples  may  be  grouped  in  any  number  of  groups.  After  the  sample  IDs 
for  the  last  group  have  been  entered,  enter  0 for  the  number  of  samples  in  the 
next  group. 
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